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Abstract 

The dynamics of the phase transition of the continuum $f+i-theory in Light Cone Quantization is 
reexamined taking into account fluctuations of the order parameter < $ > in the form of dynamical 
zero mode operators (DZMO) which appear in a natural way via the Haag expansion of the field 
$(1) of the interacting theory. The inclusion of the DZM-sector changes signiflcantly the value of 
the critical coupling, bringing it in agreement within 2% with the most recent Monte-Carlo and high 
temperature/strong coupling estimates. The critical slowing down of the DZMO governs the low 
momentum behavior of the dispersion relation through invariance of this DZMO under conformal 
transformations preserving the local light cone structure. The critical exponent 77 characterising 
the scaling behaviour at fc^ ^ comes out in agreement with the known value 0.25 of the Ising 
universality class. 77 is made of two contributions: one, analytic (75%) and another (25%) which can 
be evaluated only numerically with an estimated error of 3%. The /3-function is then found from the 
non-perturbative expression of the physical mass. It is non-analytic in the coupling constant with a 
critical exponent uj — 2. However, at D — 2, uj is not parametrisation independent with respect to 
the space of coupling constants due to this strong non-analytic behaviour. 



PACS : ll.lO.Ef, ll.lO.St, 11.30.Rd 
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1 Introduction 



Light cone quantization (LCQ) of scalar field theory in two dimensions {LC — (/>i+i) has been studied 
rather extensively during the last years. The early attempts [Q used discretization in the space direction 
thereby dealing with the problematic zero mode problem of the field operators in a simple way. Besides 
its interests as a testing ground for techniques in LCQ such scalar field theory is mainly useful for its 
phase transition properties, well known in the context of equal time quantization (EQT), and relevant to 
the description of many physical systems. 

In LCQ it is recognized now that the phase transition is driven by a coupling between two distinct 
parts of the field operator [||: the static zero mode and quantum fiuctuations of the particle field. The 
zero mode carries essential non-perturbative informations which in EQT are present in the non-trivial 
vacuum structure. It is important to realize that critical behaviour studies, were it in LCQ or EQT, can 
only be achieved for effective theories in the continuum. Hence in the discretized versions the approach to 
the continuum limit must be under control, which may be problematic |3| . Recently an LCQ- formulation 
in the continuum (CLCQ)Q was proposed in which both the zero mode and renormalization problems 
are treated in a consistant way. It is based on the treatment of fields as operator-valued distributions. 
The critical coupling gc and critical exponents were calculated to lowest order in the Haag expansion 
of the field operators. Whereas for first order perturbation theory gc comes out 45% below the value 
obtained from Monte Carlo or high temperature estimates , it is only 20% below that value for lowest 
order CLCQ. However critical exponents are usually those of mean field theory. 

To improve this situation we extend here the theory to second order in the Haag expansion of the 
fields. The benefit is a very natural inclusion of fiuctuations of the order parameter. Clearly they are 
essential to discuss the dispersion relation in the vicinity of the phase transition, a feature closely related 
to the dynamics of the phase transition and to the critical behaviour. The following results are obtained 

1) - The critical coupling is significantly improved and comes out sligthly below the canonical value. 

2) - Non-trivial values of the critical exponents are related to the existence of dynamical zero modes 
transfering momentum and not considered so far. They enter both in the equation of motion of the 
classical field (t>c{x) and of the particle field ip2{x), thereby leading to a distinctive dispersion relation. 

In Section 2 the formalism of Haag series in second order is introduced and the zero mode concept is 
extended to a dynamical operator able to transfer momentum. From the equation of motion (EM) for 
the total field ^(x) matrix elements of different types are taken giving coupled equations for the different 
components of ^{x). The solution of these equations is presented in Section 3, leading to a transcendental 
equation for the critical coupling. Its solution is discussed and solved hereafter numerically. In Section 

4 the A;^ ^ limit of the dispersion relation is discussed at some length, to extract in a transparent 
way its specific feature leading to a non-trivial value of the critical exponent 77 (ie limfc^o ^2{k) ~ k^~'^). 
Finally our results for /3-functions, critical coupling and critical exponent uj are collected and in Section 

5 a general discussion and conclusions are also given. 

2 Haag series, equations of motion 

2.1 Haag series in second order and dynamical zero modes 

The Haag series of the total field is written in terms of normal ordered products of the free field 

(po{X) with LC coordinates {X-,X+) g] 

^X)^MX) + ^{X), (2.1) 
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with 



^{X-,X+) = MX-,X+) + J dX^dX^ : MXi,X+)MX^,X+) : 

g2{X^ -X-,X^ -X-;X) + ... 
:= (po(^",^+) + 'P2(^",^+)+ higher order terms (2.2) 
in the number of free fields operators ipa. 

In Eq. (2.1) 0c(a^) and 'i'{X) are respectively the classical and quantum parts of the total field ^{X). 
Due to the interaction, 0c(-^) and ^{X) are coupled to each other. This coupling is made explicit by the 
X-dependence in the wave function (72 — X~ , X2 — X~ ; X) . It will be seen more clearly later on when 
the coupled equations for §2 and 4)c{X) are written down. It should be noticed that the contributions 
ip2[X) and (p3{X) are of different character: while (pziX) has perturbative contributions, ^P2(X) is purely 
non-perturbative and is only due to fluctuations of the order parameter. 

In CLCQ the free field ipo{X) is an operator valued distribution Q| with a generalized Fock expansion 
built out of C°°-test functions with compact support in J? ^ . It is written as 



ipoiX) = ip^{X) + ^+{X) 



dk+ 

[aik+)ei''-^ + a+ik+)e-i''-^]9ik+)fik+), (2.3) 



\/47rfc+ 



where 9{k) is the usual step function and the creation and annihilation operators a+(fc+) and a(fc+) 
obey [a{k+),a+{k'+)] = S{k+ - k'+), k.X = k-X+ + fc+X", fc" = (m^ > 0). The properties of 
the test function f{k^) are discussed in Q : its effect is to render the integral finite at both ends of the 
spectrum in and the parameter entering any possible form of f{k^) is related in an essential way to 
the final renormalization procedure of the theory. 

Written as in Eq.(2.2) ip2{X) is not apparently covariant since the integrations are done on a fixed 
time sheet. However manipulations on 0-functions, standard in Eq.(2.3), can bring it into a manifestly 
covariant form exhibiting advanced and retarded parts of the amplitude 172 ■ 

The form (2.3) is used in Eq.(2.2) and a Fourier transformation with respect to X^ — X^ and X2 —X^ 
is then performed, giving 

^,(x) = / &ieikt)9{kr)fik+)f{k+) 



"-1 "-2 

{g^+(fci,fc2;X)[a+(fc+)a+(fc+)e-^('=i+'==)-^ +a(fc+)a(fc+)e^('=i+'^=)-^] 
+ [.92+"(fci,-fc2;X)+.g2-+(-fc2,fci;X)]a+(fc+)a(fc2+)e*('=^-''^i)-^}. (2.4) 

Here ki — [kf ^ k^) with fc~ defined after Eq.(2.3). Because of hcrmiticity the coefficient fmiction 
of a(fcj*')a(fc^) is identical to 5^^. Hereafter the combinaition (7^ (fci, — ^2; X) + (7^^(— ^2, ki; X) always 
occurs in this form and will be abbreviated by G2(fci, — fc2; X). The zero mode of the field operator ip2, 
introduced in |^ where it is called Cl, is contained in the form (2.4). It is obtained after integration over 
X~ with an invariant measure which, if momenta are measured in units of the renormalized mass fi as 
in 0] (see below) is simply -^dX' . In this case, with G2 denoting the X~-independant part of G2, the 
specification is just 
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--G2{ki,-k2;X)\zM = Y+G2{ki,-kx)5{kt -k+) 



k'^ 



"-2 

:= C^{k+)5{k+ -k+). (2.5) 

It is important to realize that the zero mode concept - i.e. an X~ -independent field contribution - 
can be extended to a dynamical operator transfering momentum. The physics behind this process is as 
follows : the coefiicient function G2{ki,—k2\X) = {ki,—k2;X) + (7^~''(— ^2, fci; X) in Eq.(2.4) can 
take a periodic space dependence if the particles propagate in a periodic background field, the field (l)c{X) 
in this case. Depending upon the sign of k^ — k^ this periodic evolution is of the form exp[±|fc+X~] 
and compensates the periodic factor cxp[^(fc^ — k^)X~] coming from momentum transfer. Integrating 
Eq.(2.4) over X~ in this case generates a particular type of zero mode, which we note Cl/.- It takes the 
general form 

nk = j dk+Ci{k+,k+)a+{k+ + k+)a{k+), (2.6) 

where the coefficient Ci{k^ , k~^) is related to the coefficient function G2.o{ki , —k^ — k^) of the periodic 
part of G2{ki, —k2\ X). For fc+ = one obtains the contribution of the classical zero mode of Eq.(2.5). 
For fc"*" ^ this zero mode component of ^p2{X) has non-zero matrix elements between one particle states 
differing by a momentum transfer ±fc"'". As we shall see in this case the lowest order pertubative result 
indicates that covariance of the dispersion relation is explicit if momenta are measured in units of the 
external momentum k. A specific analysis of the infinite volume limit is then necessary, thereby leading 
to a special form of the invariant measure [rf-X""]. The specification is now 

^ G2{ku-k2;X)\DZM = , ^ g2,o(fc^ - fc+)^(fc+ - fc+ - fc+) 

V^k+{kt + k+) 

:= Ci{kt,k+)5{kt ~kt -k+), (2.7) 



\J kik2 



with V originating from the invariant measure. It is the existence of this dynamical zero mode (DZM) 
which leads to non-trivial dispersion relation and critical exponents. 

2.2 Equations of motion 

The equation of motion (EM) for the total field ^{X) is 

[49+a_ + m'] [MX] + ^{X)] + A [^^{X) + ^{X)f = 0. (2.8) 

The driving term in this equation is (j)c{X). Since we are interested in properties of the system close 
to the phase transition where all quantities related to long wave length ffuctuations of (the order 
parameter (j)o in the limit A; — > 0) are small, only terms linear in y>2 will be kept, giving a set of coupled 
equation for (f)c{X) and (/?2- 

[4d+d.+m^]MX) = -^mX)+3MXW{X) 

+ vo^{X)ip2{X) + ip2{X)ipo\X) + ipo{X)ip2{X)ipo{X)], (2.9) 
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+ MxWix) + MX)MX)MX)]. (2.10) 

Denoting the vacuum state as |0 > such that a{k'^)\0 >= 0, fc+ > 0, matrix elements of EM are 
taken between different states : < 0|£;M|0 >;< q^\EM\q^ >;< 0\EM\q^q^ > and < qlq:l\EM\0 > 
with \q^ >= a+(g^)|0 > etc.Eqs (2.9,2.10) give then a coupled system for (j)c{X),g2~^ and G2- The 
derivation is straightforward but tedious. We only give the results for the Fourier components, 4'c,gt^ 
and 6*2, with a renormalized squared-mass jjL^ 

coming from ^pQ'^{X) contributions present in the right hand side of Eqs.(2.9,2.10): 

[^,^-k+k-]Uk\k-) + ^j ^^^e{kt)e{kt)nkt)nkt) 

[\G2{kt, -k+; k) + ~gt+{kt,k+; k)\ = 0, (2.12) 

A r (7t+ 

[m' - (91 - 92 + kf\G2{qt,-q^:k) + —J --^e{k+)f{kt) 

[G2{k+, ~q+;k) + G2iqt,-kt; k) + 2~g++{kt, qt; k) + 2g++{k+,q+; k)] 

+ ^4>,{k+,k-) = 0, (2.13) 



[M^ - (91 + 92 + ^)^]^2++(9i+, 92+; ^) + / ^0{kt)f\kt) 

[G2(g+ k) + G2{qt, -kf; k) + 2g++{k+, q+; k) + 2g++{k+, q+; k)] 

+ ^4>,{k+,k-)^0. (2.14) 

It is clear that the driving term for G2 and g^"*" is (j)c{k'^ ,k~). Next a closer inspection of these 
equations shows that it is the X-dependence of 52 (cf. Eq.(2.2)) which induces fluctuations of the order 
parameter (pc. In turn this X-dcpcndcncc has its origin in the driving term (f)c present in Eqs.(2.12-2.14). 
To facilitate the analysis further it is useful to separate explicitely the zero-mode contribution of '^2{X) 
to Eqs.(2.12-2.14). However one has to distinguish between the two cases fc+ = and fc+ ^ 0. The first 
case yields the equation for critical coupling and the second the covariant dispersion relation. 
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3 Solutions of the coupled equations at = 0. Critical coupling 
3.1 Equations of motion 

At fc"*" = the explicit zero mode contributions to Eqs.(2.12-2.14) are obtained with the prescription of 
Eq.(2.5), with the notation C{g,q^) instead of C[j{g,q^). Using 

- {q, ± q,f = - {qt ± qt)i^ ± ^)] = ^ ^Jll±^^pMl , (3.I) 

9i 92 9i 92 

and with dimensionless momenta (in units of /x), Eqs.(2.13-2.14) may be reduced to (<? = 

62(91,-92) = 2ff2o(9i,-92)-§( .'T ) /^/'(fci)[2ff2++(^i,9i) 

^ 9i + 92 — 9i92 J Ki 

+ 2^2++ (fci, 92) + 62(91,-^1) +G2(A;i, -52)], (3.2) 



5^^(91,92) = g2o(9i,92) + 7( 2 , !2^1„ ^ ) / T^/^(fci)[2g^"^(fei,gi) 

4 9l + 92 + 9l92 J Ki 

+ 2g++{kuq2) + G2{qi,-ki)+G2{q2,-ki)]. (3.3) 

Here, for ease of notation, all the variables {qi,ki) stand for {qf ,k'l) and the source function g2o is 
defined as 



g2o(9i, 92) = f( ,2 , . )[/'(gOg(g, 91) + /^(92)C(g, 92) + 20o]. (3.4) 

4 9i + 92 + 9i92 



The presence of the zero mode in the source term g20 of Eqs. (3.2,3.3) is now explicit through the 
coefficients C(g, 91), C(g, 92). 

3.2 Constraints 

The projection of the EQM (2.9,2.10) onto the one particle and vacuum sectors yields the constraint 9^ 
with matrix elements < qi\03,\qi > and < 0|^3|0 >. Their vanishing give for the first one an equation to 
determine C{g, 91) in terms of g2~^ and G2 

C{g,qi)[qi+gfHqi)] + 2g J ^/2(A;i)[g2++(A:i, 9i) + ^(62(^1, -91) + 62(91, -^i))] 

+ g(po = 0, (3.5) 

and for the last one the equation (cf Eq.(2.12)) which fixes the critical coupling once §2^, G2 and C 
are known. It reads 

^o + fj ^^P{ki)nk2)[gV{k,M) + \G2{ku-k2)] + lj ^f\ki)C{g,k,) = 0. (3.6) 
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3.3 Iterative procedure 

The starting point is the analysis of Ref Q in which the terms in g^^ and G2 in Eq.(3.5) are disregarded. 
The properties of /^(gi) permit to write the lowest order solution C*^°-' as 

This fixes the initial som'cc function 520'' (qi J ^2) to be used in Eqs. (3.2-3.3). With f^{qi) = ^(91 — 
— qi), where 6{x) is the Heaviside step function, by direct computation after two iterations in Eqs. (3.2 

-3.3) it is seen that and G2 keep the same shape as 520^(91, 92) and 2g^ {qi, —q2) respectively. The 
full solution is thus seeked under the form of 3201 Eq-(3.4), with however C{g, qi) replaced by C'^^\gf, qi) 
where gf is an effective coupling to be determined. There are two reasons behind this approach: 

i) with (3.4) and (3.7) the results after one iteration for and G2 are given in closed forms. One 
finds 

hm .92 ' '{qi.qi) ^ -^gifyo ^ .9-7^]— + 0(— ), 

qi^oo b b y qi qi 

which is just 520(91,91), Eq.(3.4), in the same limit with an effective coupling gf for C^^'>{gf,qi) : gj 
= .9(1-5^)- 

Successive iterations build up the alternating sign infinite series in powers of g^^^ so that one expects 
gf to be closed to g[l + g-'^-^]^^, the departure coming from the coupling of Eqs. (3.2) and (3.3). 

ii) the full solution for C{g,qi) from Eq.(3.5) is not expected to change drastically from the lowest 
order result of Eq (3.7). Indeed using in Eq (3.5) for and G2 the ansatz build out of g2o with a not 
yet known coefficient C{g,q) one finds 

''^'^'^^-X^)-^iTq)l *^(^'^)[(fc2 + g'+M~^(fc^ + 9'-^9)^- ^'-'^ 
To a good approximation one may write 



00 



Hence with this approximation the integral term in Eq.(3.8) is exactly zero. In any case C{g,q) will 
remain close to its lowest order form (c/. Appendix A for details). 

To take into account the coupling between Eqs (3.2) and (3.3) the effective coupling gj is written 
under the form 

gf = — ^- (3.9) 

(1+5^) 

The function x{g) is expected to be close to 1 and is determined by a least-square optimisation over 
a set of grid points {qi,qj} of the solution of Eqs. (3.2) and (3.3) with the ansatzs for g^'^ and G2 as 
explained above, (cf. Appendix A). For < 17 < 7 it is found that a good fit to numerical values is 
obtained with 

^(<?) = 1 + ^. (3.10) 
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The critical coupling gc is finally obtained by evaluating Eq.(3.6). With j{k) given in term of step 
functions all the integrals are obtained in closed forms. Divergences are all proportional to ^f^{k) 
and are taken care of by a mass- type counter-term as discussed in Ref Q . The equation for is found as 

1 - |(1 + pnVs) ln[g] + j^gW3\n[l + pnV3{l + ^)] = 0. (3.11) 

To lowest order in the Haag expansion the equation was just 1 ~ f = with the solution 

ijc"^ = 4.2.... Substracting the new equation from the zeroth order one gives a small and always negative 
quantity increasing with g, indicating that > 4.2. The solution of Eq.(3.11) to second order around gc"'' 
is gc = 4.76... while the direct numerical evaluation gives gc — 4.78. For comparison with other studies 
using the convention of Parisi et al ^ this has to be translated Q| into the reduced coupling unit r which 
is just 1 at the perturbative one loop level. To the above values of gc^'^ and gc corresponds the values 
r = 1.5 and r — 1.71. This is in agreement with recent estimates from high-temperature (r — 1.754) 
[0, 1^ and strong coupling expansions (r — 1.746) ||^ and Monte Carlo simulations (r = 1.71) ||l^, p^ . 
However the recent RG-improved fifth order perturbative result of Ref. is ~ 1.837(30). A recent 
high precision estimate jl^ of 174 for the 2D Ising model leads to r = ^ = 1.7543.... The discussion on 
the peculiarity of the perturbative result is postponed to the paragraph on critical exponents. 



4 Dispersion relation, critical exponents and /3-function 
4.1 dispersion relation 

Next we discuss the dispersion relation for the order parameter fluctuations, i.e. for the field (f)c{X). Near 
the phase transition the equations of motion (2.12-14) are dominated by the dynamical zero modes. As a 
first step, in keeping with the treatment of Eq.(3.5), the lowest order solution is obtained by disregarding 
the non-zero mode contributions to Eq.(2.13). With the expressions of Eq.(2.7) for G2, Eq.(2.13) reduces 
to 



f^^G2oiqt,-qt -k+) + 



SttV 



9i 



9i 



k+ 



(4.1) 



and, with the properties of f'^{q^), 

Ci'\g,qt,k+) 



gMk)^iqt + k+)qt 
'V{q+ + k+)q+ + mq+ + k+y 



(4.2) 



Going back to Eq.(2.12) for (pdk) the DZM contribution shows up explicitely 



+ terms coming from non-zero modes = 0. (4-3) 

G20 being proportional to (f>c{k) this DZM term, irrespective of other contributions due to corrections 
arising from the coupling of with G20 through Eq.(2.14), is therefore the first non trivial correction 
to the dispersion relation of the non- interacting system (on-shell condition). 

With Eq.(4.3) as it stands the question of the covariance of the DZM contribution immediately arises. 
To examine this issue it is useful to recall how covariance emerges in the perturbative treatment with 
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LC variables. Without the zero-mode the solution for G2 and to lowest order in A is obtained from 
Eqs.(2.13-14) as 



A Mk) 
47r - (gi - 92 + m 

and 



G2o(gi, -92; k) = (4.4) 



9i iqi,q2;k) = - — --^ — ■ —rr^- (4.5) 

Sir [11^ - [qi +q2 + ky) 

Using these expressions in Eq.(2.12) the corresponding contribution is shown in Fig. 1 

1,+ _ 1,+ _ 1.+ 



k+ 


( k^ \ 











FIG. 1. Perturbative contribution to the self-energy at 2nd order in A 
With LC variables it reads 

_ f ^ r^-'" dj4 P{kt)P{kt) 

^ ' ^ ~ 96.^ h kt Jo kt - ki - kt){k- - - f - ■ 

Evaluating all momenta in units of k^ , that is k^ = xk^,k2 = yk'^ , I{k'^,k~) becomes evidently 
covariant, i.e. I{k~^,k~) = I{k^). 



[ dx [ 

967rVyo Jo [^xy{l - X - y) - xy - {x + y){l - X - y)]' 

Thus in Eq.(4-2), with q'^ expressed in units of A;+, the volume V, which has dimension of length, 
must also be written in units of -pp, i.e. V is just -p- times a function of Lorentz scalars. What is its 
argument? The external source function (jid^) (c/Eqs.(2.9-10)) actually serves to probe characteristic 
distances of the system. For a periodic background field the interaction provokes a dispersion which 
induces another length scale, the energy flow scale, which is related to the off-shellness of the process. 

It is measured in terms of k~ , say Hence in effect V oc {■^)"{^)^^~°'^ = -^{^)^^~"^ ■ Since a is 

,2 

arbitrary any scalar function of ^ is legitimate. Thus 



V 



1 e 



k- 



(4.8) 
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with v{z) — zw{l/z). The infinite volume hmit is then performed such that (k^ 



o,fc- 



oo) with 



k fixed. It is interesting to note that the two forms of V in Eq.(4.8) are related by a conformal inver- 
sion A;+ n"^ k~^ / k"^ , k~ — > ^^k~ jk"^. We shall come back to the specific LC transformation properties 
susceptible to fix further the explicit form of v{z). 

With Eq.(4.8) the DZM coeflicient cf^ of Eq.(4.2) takes now the following form, 



v{K)x{l+x) + a^{2x + l) 



(4.9) 



The zero modes covering the whole system, the presence of the volume in C^"' of Eq.(4.2) is not 
surprising. For a continuum formulation the finiteness of cf*^ in the infinite volume limit is not evident 
at first sight. In the standard discretized approach (DLCQ) all momenta are inversely proportional to 
the size L of the compact domain, indicating that cf^ should be finite in the infinite L limit. However 
the restauration of covariance in this case is not clear ||l^. It is plain to see that in the continuum ap- 
proach the infinite volume limit can be taken in a consistent way which, in keeping with the perturbative 
approach, restores also the covariance of the DZM expressions. 

Going back to Eq.(4.3) the DZM contribution can be written as lDZM{k'^)<t>c{k) 



lDZMik^)Mk) 



dx 



247r 7o ^x{l + x) 



C°Ag,x,k) 



A.9 



:ln 



The dispersion relation becomes simply 

e-fl^^lDZM{k'')^0. 

In the limit k'^ ^ it should reproduce the constraint (after renormalization) 

6 

The comparaison of these last two relations shows that limfe2^o ^i^) should be zero. For v{- 



the following expansion holds 



-^iDZAlik^ 



-^Mg)-H-) + -] + 0{v'Hv)). 



(4.10) 

(4.11) 
(4.12) 

(4.13) 



In this case the divergent part is given by the ln(w) term and is taken care of by renormalisation. If 
lim/j2^o '^(^) °^ (fc^)" with a < 1 the leading term in the dispersion relation is the one linear in v [Q. 
Apparently within our approximations for the solution of the EQM the dependence of v on fc^ determines 
completely the dispersion relation and therefore the critical exponent rj. Given that rj is sensitive to 
extensive properties of the system this result is not surprising. 

,2 

To proceed further we discuss now the physical conditions possibly fixing the function v(^). The 
volume V of Eq.(4.8) should be at least stationnary under variations of {k^,k^) preserving the local 
light cone structure. The identity map for the upper-half plane E with fc+ > is given by the projective 
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special linear group PSL{2, IR ). Let z — k + Ik^ , then an automorphism of S must be of the form 

= (4.14) 

with a,/3,7, 5 G J? and a6 — jf3 > 0. This can always be reduced by projection to a6 — = 1. 
The conformal transformation (4.14) is the product of a dilatation, a translation and an inversion. The 
translation is purely real and concerns only k^ (the energy), the origin of which is arbitrary. Since in the 
infinite volume limit dilatations do not change the physical situation special conformal transformations 
(SCT) are then the only relevant ones to consider (c/ Appendix B). SCT's corresponds to a = — 
0, 7 = c and S — 1. To derive a differential equation for v{j^) we consider the infinitesimal form of T{z) 

ST{z) = z-6cz^. (4.15) 



Putting V = V{z,z) and X ^ the evaluation of AV = V{T{z),T{z)) - V{z,z) in the limit (fc+ 
0, k~ — > oo, k^ fixed) {cf. Appendix B) , yields the differential equation: v'{X) = f^^x^ ' ^^^^ ^(^) ~ ^ 
the solution becomes v{X) = X^s) = (^)^ yielding -q = 2 — 4/3 = 2/3 for the critical exponent. For the 
universality class of the Ising model the exact result is rj = 1/4. Hence the simplest DZM contribution 
4/3 already accounts for about 76% of the exact result 7/4 . However there are corrections to the DZM 
coefficient Cf^ of Eq.(4.9) coming from interaction terms neglected in Eq.(2.13). Can they provide the 
missing contributions? With their inclusion one can see that Eq.(4.9) becomes 



C,(,, k) = - 9V^i^)[Uk)+A[g,.M 

v{^)x{l + x) + ^^{2x + l) + 2gxB{g,x,k) 

The calculation of the functions A{g^x^k) and B{g,x,k) is detailed in Appendix C. For g = 4.2 the 
coefficients Ci{g,x,k) above and cf\g,x,k) , Eq.(4.9), both divided by (l)c{k)^Jx{x + 1), are shown in 
Fig. (2. a) as a function of x for fc = 0.1/i and in Fig. {2.h) as a function of k for x — \. 




FIG. 2. Ci{g,x,k) and C} {g,x,k) as a function of x for k = 0.1 (2. a) and of k 
for a; 1 (2.5), both for g = 4.2. 
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With this new value of Ci the dispersion relation integral lozAiik'^) can only be evaluated munerically 
for given values of g and k. Fig. (2. a) indicates that the result should be close to the initial one with 
however some difference in the fc-dependence as seen from Fig. (2.5). In the limit fc ^ and for g = 4.2 
the initial dispersion integral lozuik'^) = ^^lozAiik'^) of Eq.(4.10) reduces to (c/ Eq.(4.13)) 

lim iozMik^) = -2.12585 + 1.3333 log(fc) - 0.238663/c(3) + 0.046304A:(tl (4.17) 



fc^o 



The two functions lozMik"^), evaluated as above and by numerical integration with Ci{g,x,k), are 
plotted in Fig (3. a) as a function of k iov g = 4.2. In the later case a least square fit (c/. Appendix C) 
with a functionnal of the form of Eq.(4.17) gives 

VrniiDZMik"^) =-2.3278+ 1.3333 log(fc)- 27.103fc(5)+ 3040.6fc(5\ (4.18) 



Ci 



which is compared to the result of the numerical integration with C\ in Fig. (3.b). 















CIO = 




/ 




CI = 



















£lt = 






exact - 





(3.a) 



(3.6) 



FIG. 3. The dispersion integrals lozMik^) 

numerical and fitted {Eq. {4.18)) expressions in (3.6) as a function of fc for 5' = 4.2. 



and lozMik"^) in (3. a) and the 

C° Ci 



It seems remarkable; that the expected power of fc, | giving 77 = 2 — | = i, provides quite a good 
fit. However, despite the confidence in our fitting procedure -it has been tested in fitting the exact 
integration performed with c[^^ instead of Ci, whereby the result (4.17) is reproduced with high accuracy 
(c/. Appendic C) - the least square criterion is not stringent enough to fix the power of k with total 
trust. E.g. a 10 % variation of the per point results from the arbitrariness of the distribution of points 
between k=0 and k=0.1. This leads to a corresponding uncertainty of the power of k of about 3% which 
locates 7] between 0.2 and 0.3 or in other terms 77 = 0.25 ±0.05. An analytical analysis of the behaviour of 
the full dispersion relation at low fc^ is the only way to settle the issue without ambiguities. Unfortunately 
we have not been able to perform such an analytical investigation of the interaction corrections to the 
dispersion integral lozAtik'^)- 



4.2 /3-function 

To determine the /3-function the constraint, Eq.(3.5), which depends on the cutt-off A through the test- 
function f^{k), is considered as a prescription for the calculation of the critical mass M{g,A). M'^{g,A) 
is just /U^ times the right hand side of Eq.(3.11), with /x^ = X/{A'Kg). 
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One has 



dg 

Evaluating (4.19) gives 



/3(ff)=M[— ]^a'a)- (4-19) 



with, 



Pig) = -2.9^> (4.20) 



N r, 9/-, 27rv3, , , . ,7rv3, 7rv3/, 9 xmt, ""VO/, 9 

^(5) = [l-|(l+5-^)ln(5)+3'^ln(l+5^(l + ^))][l+5-^(l + ^)], 

i5(g) = (i + 5)[i + ^W3(i + A)] + 5W3 gW3 
u) V -r gn g V -r -^2''J ^ 27 ^ 108 ^ 



From standard rcnormahzation group analysis one must have 

\im dig) = -{A ^D)g + 0{g^) 

which is indeed the case for Pig) of Eq. (4.20). To lowest order in the Haag expansion the corresponding 
/3-function is just 

M9) = -^^:^. (4.21) 
(1 + fi) 



The two functions Pig) and Poig) are plotted in Fig. 4. 




FIG. 4. The functions ^(5), Eqs. (4.20-4.21). 

In both cases the exponent w = P'igcr) — 2 exactly In particular the authors of Ref |pl| emphasize 
that dX D = 2 perturbative estimates of the critical value of g and of oj are not reliable due to strong 
nonanalytic corrections to the /3-function. Moreover these nonanalyticities invalidate the conventionnal 
wisdom that oj is parametrization indepcndant with respect to the space of coupling constants. This is 
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most easily seen for the Ising model. With b denoting the usual dimensionless Ising coupling, the critical 
value is be = ^ log(\/2 + 1). Let i = 1 — then the "mass" , i.e. the inverse ol the correlation length ^, 
is a known lunction oi t, m ~ t. Using Eq.(4.20) gives = —bet and u; — (3'{bc) = 1. A field theoretic 
definition of the coupling jl^ is now introduced 

5* = lim5(t) = lunh^], (4.22) 

Here x is the "magnetic" susceptibility, is the 4-point function at zero momentum and ^ the 2nd 
moment of the correlation length. With this definition the new beta function becomes 



which implies 



/3'(g*)=/3'(6)+/3(6)|l||. (4.24) 



The common wisdom is then that a) — I3'{g*) = P'ibc) = ^ since /3{bc) = 0. However it is not so since 
g*{b) = with A 7^ 1. Indeed g*" {b)/g*'{b) = -(A - l)/(6ci) and since /?(&) = -b^t the correct relation 
is now 

Cj^uj + A-1 = A. (4.25) 
For the Ising model A — j, and for with g = uj = 2. 

The situation was apparently known long ago |l9| and is examplified here in a simple and transparent way. 



5 Conclusion 

In the framework of LC quantization the mechanisms of symmetry breaking have attracted a lot of interest 
because of their evident nonperturbative aspects. In this respect the analysis of critical behaviour is 
a key issue to which, however, little attention has been given in the past. In particular for D — 2 
it is well known that critical exponents are truely nonperturbative extensive quantities which are not 
easily calculable by conventional methods. 2D Conformal Field Theory (CFT) was the most spectacular 
achievement providing a wealth of analytical results for a variety of models. Some specific deformations 
of certain 2D CFT give a hint to the treatment of some particular 3D systems but the systematics is, so 
far, elusive po| . 

The LC quantization scheme is an interesting nonperturbative alternative valid in any dimension. 
In this paper we have shown that a full understanding of the physics of zero modes combined with 
non-compact continuum field dynamics is necessary to correctly evaluate the critical coupling and to 
treat covariantly the low k"*"— region determining the scaling behaviour. It is found that conformal 
transformations preserving the local light cone structure govern the scaling behaviour of the two-point 
function and provide 75% of the total contribution to the critical exponent rj . The remaining contributions 
being due to specific interaction terms were analyzed and found to have the expected magnitude to make 
up for the missing 25%. However, our result for these corrections is based on a numerical analysis 
and should be supported in addition by analytical investigations taking full advantage of the conformal 
symmetries valid at the critical point. Thereby one could possibly relate the CLCQ approach to the 
corresponding CFT. It is worth noting that the ensemble of multicritical scalar models (the so called 
unitary A- A series in CFT language) with effective Lagrangians £ = ^{d(p)^ + ip'^'^, k > 2, can also 
be treated in CLCQ. For k > 2 the equations of motion, the constraints and the dispersion integral 
lDZM{k^) will differ from the ones for if'^. Comparing the critical exponents rj with exact CFT results 
would provide another test of the CLCQ DZM treatment. 
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Appendix 



A Approximate solution of the coupled equations (3.2-3.3) 



For the sake of clarity a more elaborate notation is used here. The zero order coefRcient Cq of Eq.(3.7) 
depends also on the coupling g and on the order parameter so it is noted now 



Co{g,(j)o,k) = - 



{g + ky 



(A.1) 



To make clear that the source term g2o in Eqs. (3.3-3.4) may depend on the initial and effective coupling 
g and gf through the coefficient Co, 520 is noted 



52o(5,5/,</'o,«i,g2) 



91 92 



4 {qt +qI + Q1Q2) 



[Co(5/,'/'o,«i) +Co{gf,(j)o,q2) + 2(j)o]. 



(A.2) 



Then 321(3, .g, (/'o, 91) 92) and 021(5,5, ?^o, 91,92) will be the one time iterated value of 5^''' and 62 • They 
write explicitely 



521(5,5,^0,91,92) = 520(5,5,^0,91,92) - 5-; 



91 92 



520(5,5, '/'cfei, 91 ) 
(9i + 5i + 9192) L"' Jo (fc?+9?-fci9i) 



91 / dki 



J 

Jo 



+ 92 / dki 



92oig,g,(l>o,h,q2) _ 1 

{kf +ql- kiq2) 2 



r°° dki 



520(5, 5, 9^0, fci, 91) 



^ rdki ■ 

2 J -^520(5, 5, 00, Ki, 92) 



(A.3) 



91 92 



^21(5,^0,91,92) = 2ff2o(5,'Ao,9i,-92) + 25 2 , 2 ^ 

(,9i + 92 - 9i92j 



91 



i 



(^fc^.^2o(5,5,</'o,^i,9i) 



r 520(5,5,00,^1,92) 

[kf + 92 - ^192) 2 



1 p dfc 
91 



{kf + qf- km) 

-^520(5,5,00, fcl, 91) 



1 



92 

-^520(5, 5, 00, fci, 92) 



(A.4) 



The integration limits are dictated by the fact that only states with positive momenta can be created. 
The integrals in Eqs.(A.3-A.4) can be obtained in closed forms. To show that 52o(5,5/,0o,9i,±92) 



with 5/ = 



(i+s^) 



leads to good approximate solutions to the coupled integrals equations (3.2-3.3), 



522(5, 5/, 00, 91,92) and 622(5, 5/, 0o, 91,92) will denote the right hand side of Eqs. (A. 3) and (A.4) re- 
spectively with 320 in the integrals replaced by 32o(5,5/,0o,9i,9i)- The functions 520(5, 5/, 0o, 9i, 9i), 
521(5, 5, 00, 9i, 9i) and 522(5, 5/, 0o, 9i, 9i) are plotted in Fig.(A-l) as function of qi for ^0 = 0.01 and 
5 = 1.5 and 5 = 4.2. 
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g-1.5 



g=4.2 





FIG. A-1. Approximate solutions of Eq.(3.3) as explained in text. 

In Fig.(A-2) the functions 2g2o{g,gf,(f>o,qi,-^-0), G21 (5> 5> </'o, 91, 1-0) and 6*22(5, 5/, </'o, 9i, 1-0) are 
represented as function of qi for the same values of (j)o and g (purely diagonnal cases are not relevant 
here) . 



g=i-5 



9=4.2 



-0.006 - 





FIG. A-2. Approximate solutions of Eq.(3.2) as explained in text. 



In both cases the one time iterated functions keep the same shape as the source term. Other non- 
diagonnal cases in show similar features. The agreement between 520(5,5/, (Ao, 9i, 92) and 522(5,5/, (t>o,Qi,Q2) 
is already quite satisfactory and can be further improved for 2g2o(5, 5/, "Ao, 9i, ~92) and (^22(5, 5/, <p0y 9i, Q2) 
in the following way. The effective coupling 5/ is changed to 



9 fix) = 



(1+5^) 



(A.5) 



Then a set of grid points in q is chosen on which a x^-function depending on {g, x, (po) is built out of 
the difference squared between the ansatz solutions based on ^20 ,522 and G22- Thus, for a given {g, ({>q), 
X is the value which minimizes in a least square sense the difference between the left and right sides of 
the integral equations with 5^'^ and G22 replaced by the ansatz just mentionned. The following table of 
values is obtained {(po = 0.01, the precise value is irrelevant) 



g 


.2 


1 


2 


3 


4 


5 


6 


7 


x(g) 


.822 


.918 


1.0 


1.09 


1.17 


1.255 


1.33 


1.41 
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As expected X remains close to one with a linear dependence 0.0835(7 ~ For 9 small (< 2 ) the 
values of x depend somewhat on the repartition of the grid points and the analytic argument of section 
(3.3) indicates that a; = 1 in this cases. For g < 7 the form used hereafter is 

x{g) = l + ^. (A.6) 

The quality of the ansatz solutions can be vievwed easily by plotting g2o{g, gf{x), (f>o,Qi, Q2), 
522(5, <?/(a;), 00: 917 92) etc.. since all integrals in Eqs.(A.3-A.4) can be performed. 

The expressions for ^^"^ and G22 are now used to find the new C{g, <f)o,k), solution of Eq.(3.5). The 
evaluation is simple and consists in performing exactly the integrals in Eq.(3.8). For g <1 the plots do 
not distinguish between Co{g, (t>f), k) of Eq.(A.l) and the new function of k. Their difference increases 
with g but remains small as shown in Fig. (A. 3) for g = 5. 



— 0.003 -t 




2 4 , e s 

l<~ 



FIG. A.3. the functions C{k) of Eq.(3.7) and solution of Eq.(3.5) for g = 5 

This establishes the consistency of the procedure so that no further iteration seems needed with a 
more accurate expression for C{g, (t)o,k). 

B Behaviour of the volume under transformations of the pro- 
jective special linear group PSL(2, IR ) preserving the local 
light cone structure 

Let z = k~ + Ik~^. The conformal transformation T{z) = with real coefficients a,/3, 7, 6 such that 

a6 — l3'y = 1, maps the upper-half plane onto itself and preserves the local light cone structure, namely 

{timelikevector — > timelikevector 
nullvector — > nullvector 
spacelikevector — > spacelikevector 
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In full generality one may choose a = e^(l + ac), (3 — ae^^ , 7 = ce^ and S — e^^. T(z) results then 
from the successive actions of a dilatation, z ze^, followed by a special conformal transformation 
(SCT), z — » j^^t and a translation, z — > z + a. For infinitesimal values of the parameters Sa, 6c and 6p 

T{z) = z + 5a — Scz"^ + SpZ. Requiring stationnarity of the volume V = prt'(^) = V{z, z) under the 
infinitesimal transformation z T{z) leads to the relation 

+ ^-^^l^ - ^Jk^ - V + ^^^^ - Jk^^^-djT' ^^-^^ 

where X = j^. In the limit — > 0, fc~ — > 00) with fixed only the 5c (SCT) contribution survives. 
With the initial condition v{0) = the solution is v{X) = X^i\ 

C Corrections to the dispersion relation and power law be- 
haviour at small k^. 



The starting point to evaluate corrections to c[^\g, q'^ , k~^) of Eq.(4.2) is Eq.(2.13). Using Eq.(2.7) to 
extract the DZM contribution, the complete equation for Ci(g, q"*", fc"*") now reads 



+ 



Ci(g,a:, fc) 
Vx{x + l) 

^ f '-P{z) [g^^(z, -(1 + x); k) + G^^ix, -z; k) + 2g++{z, x; k) + 2~g++{z, 1 + x; k) 



)x(l+x) + ^(2x + l) +gMk) 



(C.l) 

= 0, 



where = xfc"*" and is the non-zero mode contribution to G2. The corrections come from the 

integral terms which are evaluated with the lowest order solution g^^^ and 620'^, i.e. obtained by 
neglecting integral terms in their own equation. These solutions read 



920~ix,y;k) = h{x,y]k) J-^Ci{g,x,k) + J--^c'f\g,y,k)+2(l)c{k) 

LVl + x Vl + y 



+ y 



G20 {x,y;k) = 2(j)c{k)h{x,y;k), h{x,y;k) = - — 



xy 



(C 2) 

(1 + A;2x)2/(2/ + x(l + x)) + x(l +x) ■ ^ 

In g2o'{x, y; k) terms which can be regrouped later on with C\{g, x, k) in Eq.(C.l) have been singled out 
whereas Ci ((/, y, k) was replaced by C^°^ [g, y, k) whenever this coefiicient enters under the integral to be 
performed. Eq.(C.l) can then be written as 



Ci(g,x, fc) 
^/x{x + l) 



)x(l + x) + |(2x + 1) + 2gxB{g, x, fc)J + 5(<^c(fc) + A{g, x, k)) = 0, 



(C.3) 



where A{g, x, k) = T^f^-^^Ai^g, x, k) + 24)c{k)[B{g, x, k) + B{g, 1 + x, k)] with 



Ai{g,x,k) 
A2{g,x,k) 
As{g,x,k) 
A4{g,x,k) 



= \r ^G2T(.,-(l + x);fc) 

^ -'(l+a;) ^ 

— o I — G2o^ {x,—z;k) 

= / —\T^'^'^°''(9,z,k)h{z,x,k) 
Jo z y L + z 

f°° dz 

= A3{g,l + x,k), B{g,x,k)= / — 

^0 z 



h{z, X, k). 



(C.4) 
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All the integrals in Ai{g, x, k) and in B{g, x, k) can be obtained in closed forms but the final integral over 
X giving the dispersion contribution lozMik"^) can only be done numerically for given values of g and k. 
As seen in Fig. (2. a) Ci{g,x,k) and c[^\g,x,k) approach each other for large x; we therefore separate 
the integral over x as follows 

/ / . . '^i{9,x,k) = [ Ci{g,x,k) + 

Jo ^/x{l + X) Jo ^yx(l+X) 

CMk) r- dx ^(0) 

(5,50,fc)i50 ^/x{l + x) 1 ' ^ ^ > 

Thereby we avoid delicate numerical work in the limit of large x, since the last contribution can be 
evaluated analytically. For g = 4.2 the corresponding dispersion integral lozMik^) is shown by the 
continuous curve in Fig. (3). A fitting procedure is now envisaged to extract its power law behaviour as 
a function of k in the small k region. The method is tested first on the known result with c[^\g,x,k) 
which reads 

VirnlDZMik^) = -2.12585 + ^ log(A:) - 0.238663A;(3) + 0.046304A:(t) - 0.037973fc(t) log(fc). (C.6) 

A grid of 40 points between fc = and 0.1 is chosen according to the relation 

k{i) = i tan[i arctan[0.210](l + z{i))], (C.7) 

where z{i) arc the 40 Gauss points with — 1 < < 1. This choice produces points with greater density 
in the region where the function to be fitted is varying most rapidly. With a trial function of the form of 
Eq.(C.6) with unknown coefficients and powers of k the constant and the coefficient of the log term are 
first adjusted by least-square fitting the lower part of the spectrum. Thereby the domain to explore in 
the multiparameter space is considerably reduced. A full least-square search is then performed over the 
whole 40 points giving 

fit{k) = -2.12585+ 1.33333 log(/c) - 0.23848fc(^-^-''^''-'') + 0.05382fc(=^-^*'^*'' - 0.03335/c(^-^^'^®' log(fc). (C.8) 
The same tactic is used for the numerical spectrum of Fig. (3). The spectral values are given in table (C.l) 



k 








k 






k 




i(k) 




k 






i(k) 


9.12 X 10" 


5 


-5.484 


5.396 X 10 


-2 


-3.36 


5.396 X 10 


-2 


-3.36 


9 


035 X 


10-^ 


2 


-2.855 


4.799 X 


10" 


-4 


-4.816 


5.801 X 10 


-2 


-3.317 


5.801 X 10 


-2 


-3.317 


9 


308 X 


10" 


2 


-2.802 


1.177 X 


10" 


-3 


-4.472 


6.204 X 10 


-2 


-3.272 


6.204 X 10 


~2 


-3.272 


9 


556 X 


10" 


2 


-2.752 


2.178 X 


10- 


-3 


-4.252 


6.6 X 10" 


2 


-3.224 


6.6 X 10" 


2 


-3.224 


9 


778 X 


10" 


2 


-2.705 


3.477 X 


10- 


-3 


-4.098 


6.99 X 10" 


2 


-3.175 


6.99 X 10" 


-2 


-3.175 


9 


972 X 


10" 


2 


-2.661 


5.066 X 


10- 


-3 


-3.986 


7.369 X 10 


-2 


-3.124 


7.369 X 10 


-2 


-3.124 


1 


014 X 


10" 


1 


-2.623 


6.937 X 


10- 


-3 


-3.9 


7.735 X 10 


-2 


-3.071 


7.735 X 10 


-2 


-3.071 


1 


027 X 


10" 


1 


-2.59 


9.077 X 


10- 


-3 


-3.833 


8.088 X 10 


-2 


-3.017 


8.088 X 10 


-2 


-3.017 


1 


038 X 


10" 


1 


-2.563 


1.147 X 


10- 


-2 


-3.78 


8.423 X 10 


-2 


-2.963 


8.423 X 10 


-2 


-2.963 


1 


045 X 


10" 


1 


-2.545 


1.411 X 


10" 


-2 


-3.735 


8.739 X 10 


-2 


-2.909 


8.739 X 10 


-2 


-2.909 


1 


049 X 


10" 


1 


-2.534 



Table C.l. Values oilozMik'^) 



as function of k 



21 



A least-square fit on the eight first points of the spectrum fixes the constant and the coefficient of the 



log term to —2.185 and 1/3 respectively. The remaining part lozMik'^] 
Fig.(C.l) and is therefore a power of k close to 2. 



Ci 



+ 2.185 — hlog{k) is shown in 




FIG. C.l. iDZAiik^) + 2.185 - ^log{k) as a function of k. 

Ci 

The final least-square search with a trial fitting function a + ^log{k) + ck^ + dk'^'^ is made over the 
whole spectrum with input values a = —2, c = 27, e = 1 and d = 3000 with the result 

fit{k) = -2.3278 + \ log(fc) - 26.9999fc(i-^^9^9) + 3040k^^-^^^^l (C.9) 

Fixing finally e to | gives the fit of Eq.(4.18). 
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